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higher  than  1500  m  are  shaded. 
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events  described  in  Table  2. 

Paths  used  for  single-station  studies  of  multi-mode  Rayleigh-  8 
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(Top)  Original  trace  (in  counts)  for  the  earthquake  that  occurred  12 
on  January  28,  1996  as  recorded  at  station  SODA.  (Bottom)  Win¬ 
dowed  trace  of  groimd  displacement  sifter  correction  for  geometri¬ 
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Contoured  spectral  amplitudes  of  fundamental-  and  higher-mode  13 
Rayleigh  waves  obtained  for  an  earthquake  that  occurred  on  22 
April  1995  and  recorded  at  station  GNI.  The  windowed  ground 
displacement,  after  correction  for  the  geometrical  spreading  and 
instmment  response,  and  bandpass  filtering,  is  displayed  at  the 
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Interstation  group  (bottom)  and  phase  (top)  velocities  for  the  sta-  15 
tion  pairs  in  Region  1.  The  caption  shows  symbols  used  for  vari¬ 
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Abstract 


Observed  velocities  and  attenuation  of  fundamental-mode  Rayleigh  waves 
in  the  period  range  7-82  sec  were  inverted  for  shear-wave  velocity  and  shear- 
wave  Q  structure  in  the  Middle  East  using  a  two-station  method.  Additional 
information  on  Q  structure  variation  within  each  region  was  obtained  by  study¬ 
ing  amplitude  spectra  of  fundamental-mode  and  higher-mode  Rayleigh  waves.  We 
obtained  models  for  the  Turkish  and  Iranian  Plateaus  (Region  1),  areas  sur¬ 
rounding  and  including  the  Black  and  Caspian  Seas  (Region  2),  and  the  Arabian 
Peninsula  (Region  3).  The  effect  of  continent-ocean  boundaries  and  mixed  paths 
in  Region  2  may  lead  to  unrealistic  features  in  the  models  obtained  there.  At 
lower  crustal  and  upper-mantle  depths,  shear  velocities  are  similar  in  all  three 
regions.  Shear  velocities  vary  significantly  in  the  uppermost  10  km  of  the  crust, 
being  3.21,  2.85,  and  3.39  km/s  for  Regions  1,  2,  and  3,  respectively.  Q  models 
obtained  from  an  inversion  of  interstation  attenuation  data  show  that  crustal 
shear-wave  Q  is  highest  in  Region  3  and  lowest  in  Region  1.  Q  s  for  the  upper  10 
km  of  the  crust  are  63,  71,  and  201  for  Regions  1,  2,  and  3,  respectively.  Crustal 
Q's  at  30  km  depth  for  the  three  regions  are  about  51,  71,  and  134.  The  lower 
crustal  Q  values  contrast  sharply  with  results  from  stable  continental  regions 
where  shear-wave  Q  may  be  as  high  as  a  thousand  or  more.  These  low  values 
may  indicate  that  fluids  reside  in  faults,  cracks,  and  permeable  rock  at  lower 
crustal,  as  well  as  upper  crustal  depths  due  to  convergence  and  intense  deforma¬ 
tion  at  all  depths  in  the  Middle  Eastern  crust. 
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1.  Introduction 


Seismic  surface  waves  have  been  widely  used  to  study  the  shear-wave  veloc¬ 
ity  distribution  in  the  Earth.  In  recently  years,  high-quality  data  and  increasing 
numbers  of  favorably  located,  modem,  digital  stations  have  enabled  seismologists 
to  apply  tomographic  techniques  to  surface  dispersion  data.  Numerous  such  stud¬ 
ies  have  provided  information  of  regional  shear-velocity  variations  at  both  global 
and  regional  scales  in  the  Earth. 

Studies  of  surface-wave  attenuation  have  been  less  numerous  mainly 
because  of  well-known  difficulties  in  measurement;  surface  waves  may  be  con¬ 
taminated  by  lateral  refraction  and  multipathing  that  may  focus  or  defocus 
energy  and  produce  spurious  results.  It  is,  however,  important  to  study  attenua¬ 
tion  because  it  can  provide  information  on  temperature,  fluid  content,  phase 
change,  and  density  of  solid-state  defects  in  the  crust  and  mantle,  phenomena 
that  are  not  easily  studied  using  only  seismic  wave  velocities. 

Few  models  of  shear-wave  Q  (Q^)  in  the  crust  and  upper  mantle  have  been 
obtained  for  continental  regions.  But  these  have  given  us  some  information  on 
low  related  to  the  tectonic  history  of  various  regions  (Mitchell,  1995).  That 
work  suggested  that  the  value  of  obtained  in  the  upper  cmst  of  any  region  is 
directly  proportional  to  the  time  that  has  elapsed  since  the  most  recent  episode  of 
major  tectonic  activity  there. 

The  present  study  measures  Rayleigh-wave  dispersion  and  attenuation  in 
the  Middle  East,  a  complex  region  of  interaction  between  the  Eurasian,  Arabian, 
and  African  plates.  We  will  invert  two-station  data  for  models  of  shear  velocity 
(A)  and  for  three  broad  regions  and  study  smaller-scale  variations  of 
within  each  region  using  a  single-station  method.  The  models  will  then  be 
compared  to  the  geology  and  tectonic  history  of  different  portions  of  the  Middle 
East. 
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2.  Overview  of  Middle  East  Tectonics 

The  Middle  East  is  a  tectonically  complex  region  where  the  Arabian, 
Afiican,  and  Eurasian  plates  meet.  Figure  1  shows  the  relative  movement  of  the 
plates,  the  major  tectonic  features,  and  the  boimdaries  of  tectonic  provinces  in 
that  region. 

The  Arabian  plate,  consisting  of  a  shield  in  the  west,  a  central  shelf,  and 
platform  and  basins  in  the  east,  is  the  most  stable  portion  in  this  region.  It  is 
bounded  on  the  northwest  by  the  Dead  Sea  and  East  Anatolian  faults,  and  on  the 
northeast  by  the  Zagros  thrust  belt.  The  North  Anatolian  fault  extends  in  an 
east-west  direction  across  almost  all  of  Turkey.  The  Arabian  shield  consists  pre¬ 
dominantly  of  Precambrian  gneiss  and  metamorphosed  sedimentary  and  volcanic 
rocks  that  have  been  intruded  by  granites  (Powers  et  al.,  1966;  Brown,  1972). 
The  thin  basaltic  lava  flows  that  have  been  emplaced  along  the  western  border  of 
the  shield  since  Miocene  time  reflect  the  intensity  of  volcanism  in  the  region 
(Brown,  1972;  Coleman,  1977,  Mokhtar  and  Al-Saeed,  1994).  East  of  the  shield, 
the  shelf  is  tectonically  stable.  A  cross-section  through  the  peninsula  (Seber 
et  al.,  1997),  trending  roughly  east-west,  shows  that  the  thickness  of  the  sedi¬ 
ments  increases  gradually  across  the  shelf  area,  more  rapidly  to  the  east,  reach¬ 
ing  about  10  km  in  the  basin  region. 

The  Arabian  plate  and  Eurasian  plate  are  currently  colliding  along  the 
Zagros  suture  zone  in  western  Iran  (Stocklin,  1974;  Sengbr  and  Kidd,  1979; 
Stoneley,  1981).  The  Turkish  and  Iranian  Plateaus,  resulting  from  the  collision 
between  the  Arabian  and  Eurasian  plates,  are  about  1500  km  high  in  average 
(Sengor  and  Kidd,  1979).  After  continental  collision  during  the  late  Miocene,  con¬ 
vergence  between  Arabia  and  Eurasia  is  continuing  as  shown  by  the  Pliocene  to 
recent  folding  and  thrusting  in  the  Zagros,  the  border  folds  of  southeastern 
Turkey  (Ketin,  1966),  and  the  present  diffuse  seismicity  of  the  entire  Turkish  and 
Iranian  Plateaus  (Canitez  and  Ucer,  1967;  McKenzie,  1972;  Nowroozi,  1972; 
Berberian,  1976). 


3 


3.  Previous  Seismological  Work  in  the  Middle  East 

Although  the  tectonics  of  the  Middle  East  is  well  understood,  only  a  few 
reliable  models  of  shear  velocity  and  attenuation  have  been  obtained  there.  In  an 
early  study,  Niazi  (1968)  inverted  short-period  surface-wave  group  and  phase 
velocities,  and  obtained  shear  velocities  of  3.4  km/s  for  the  upper  crust,  3.6  km/s 
for  the  lower  crust  and  4.6-4.9  km/s  for  the  upper  mantle  in  the  Arabian  Penin¬ 
sula.  Knopoff  and  Fouda  (1975)  inverted  the  phase  velocities  of  Rayleigh  waves 
for  the  period  range  22-167  sec  to  obtain  shear  velocity  models  using  a  two-sta¬ 
tion  method.  The  paths  connecting  station  pairs  (SHI-HLW  and  SHI-JER)  are 
both  located  in  the  northern  part  of  the  peninsula,  trending  east-west.  Their 
study  showed  a  variety  of  solutions  which  can  explain  the  data.  Mokhtar  2md  Al- 
Saeed  (1994)  inverted  group  velocities  of  fundamental-  and  higher-mode  Rayleigh 
and  Love  waves  for  velocity  structure  in  the  Arabian  Peninsula.  They  studied 
paths  between  many  epicenters  and  station  RYD,  located  in  the  center  of  Arabian 
Peninsula.  They  divided  their  data  into  three  major  groups  and  obtained  differ¬ 
ent  models  for  each  group.  Ghalib  (1992)  used  data  recorded  at  stations  in  the 
northern  part  of  the  Arabian  Peninsula  to  invert  for  the  velocity  structure.  The 
only  models  of  velocity  structure  known  for  the  Turkish  Plateau  are  those 
obtained  by  Mindevalli  and  Mitchell  (1989)  using  the  single-station  measure¬ 
ments  of  Rayleigh  and  Love  wave  group  velocities  in  the  period  range  8-50  sec. 
Different  models  were  obtained  for  eastern  and  western  Turkey.  In  this  study, 
detailed  velocity  models  will  be  provided  for  different  regions  of  the  Middle  East. 

Q  studies  in  the  Middle  East  have  indicated  that  both  regional  phases  and 
surface  waves  attenuate  rapidly  with  distance.  Kadinsky-Cade  et  al.  (1981)  found 
that  regional  phases  were  extinguished  along  some  paths,  but  not  along  others. 
Numerical  values  for  Lg  Q  in  Iran  were  obtained  by  Nuttli  (1980)  and  by  Jih  and 
Lynnes  (1992).  Both  studies  indicate  that  Lg  Q  is  low  in  that  region,  and  the  lat¬ 
ter  study  showed  that  it  is  also  quite  variable  for  different  paths  across  the  Ira¬ 
nian  Plateau.  Seber  and  Mitchell  (1992)  inverted  surface  wave  attenuation  data 
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to  obtain  the  models  for  the  upper  crust  of  the  Arabian  Peninsula.  They  found 
that  shear-wave  Q  values  in  the  crust  vary  from  about  60  along  the  margin  of 
Red  Sea  to  100-150  in  the  central  part  of  the  peninsula,  and  65-80  in  the  eastern 
folded  region.  Ghalib  (1992)  applied  the  extended  coda-Q  technique  to  short- 
period  Lg  coda  waves,  obtaining  a  imiform  distribution  of  coda  Q  with  an 
average  value  of  214  for  the  Arabian  plate. 

4.  Regionalization  and  Data  Acquisition 

In  this  study,  we  used  both  two-  and  single-station  methods  for  both  velocity 
and  attenuation  studies.  In  order  to  model  the  regional  variation  of  shear  veloci¬ 
ties  and  attenuation  in  the  Middle  East,  we  have  divided  it  into  three  regions, 
based  upon  the  availability  of  data  for  two-station  studies.  Region  1  consists  of 
the  Turkish  and  Iranian  Plateaus  and  includes  the  station  pairs  GNI-KTV  and 
GNI-ANTO  (Figure  2).  Those  paths  are  confined  strictly  to  continents.  Paths 
crossing  the  Black  and  Caspian  Seas  are  excluded  from  this  region  and  are 
assigned  to  Region  2.  The  station  pairs  for  that  region  are  ANTO-KIV,  GNI- 
ABKT,  and  KIV-ABKT  that  cover  the  two  seas  plus  a  large  land  area  around 
them.  Region  3  is  the  Arabian  Peninsula.  For  that  region  we  used  data  from  sta¬ 
tions  within  Saudi  Arabia  as  well  as  one  station  in  Israel.  The  station  pairs  for 
Region  3  are  BGIO-RAYN,  BGIO-HALM,  AFIF-SODA,  HALM-RANI,  and  HALM- 
SODA. 

To  improve  regional  coverage  we  used  a  multi-mode  method  that  requires 
only  a  single  station  (Cheng  and  Mitchell,  1981).  Figure  3  shows  the  additional 
paths  and  improved  regional  coverage  made  possible  by  the  single-station 
method. 

The  stations  used  in  this  study  include  four  Incorporated  Research  Institu¬ 
tion  for  Seismology  (IRIS)  stations  (ANTO,  ABKT,  GNI,  KIV),  one  Geophone  sta¬ 
tion  (BGIO),  and  five  stations  in  Saudi  Arabia  (RANI,  RAYN,  HALM,  AFIF, 
SODA),  deployed  temporarily  by  the  University  of  California  at  San  Diego  and 
Boise  State  University.  The  locations  of  the  stations  are  given  in  Table  1.  All  of 
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Table  1.  Station  Information 


Station 

Latitude 

(°N) 

Longitude 

(0  E) 

ABKT 

37.9304 

58.1189 

ANTO 

39.8689 

32.7936 

BGIO 

31.7220 

35.0878 

GNI 

40.0530 

44.7240 

mv 

43.9562 

42.6888 

AFIF 

23.9310 

43.0400 

HALM 

22.8454 

44.3173 

RANI 

21.3116 

42.7761 

SODA 

18.2921 

42.3769 

RAYN 

23.5220 

45.5008 

the  stations  provided  broadband  digital  data. 

In  the  two-station  method,  27  earthquakes  were  used  as  indicated  in  Table 
2.  The  dates,  origin  times  and  locations  of  the  earthquakes  and  recording  sta¬ 
tions  are  provided  by  the  IRIS  DMC.  All  but  two  of  the  earthquakes  are  less  than 
50  km  in  depth.  The  two  exceptions  (events  20  and  24  in  Table  2)  have  depths  of 
75  and  62  km.  In  the  study,  all  events  lay  within  8°  in  a2dmuth  from  great-circle 
paths  between  the  two  recording  stations.  The  first  station  of  each  station  pair  in 
Table  2  is  the  closer  of  each  pair  to  the  epicenter.  The  rightmost  coliunn  in  Table 
2  indicates  the  region  that  the  surface  waves  traverse.  Earthquakes  26  and  27 
are  not  used  for  determinations  of  properties  within  any  region  since  the  waves 
from  those  events  traverse  three  regions.  They  are  only  used  to  see  if  velocities 
and  attenuations  along  those  paths  are  consistent  with  our  results  for  the  indi¬ 
vidual  regions. 

18  earthquakes,  listed  in  Table  3,  were  analyzed  using  the  single-station 
method.  All  but  one  of  the  earthquakes  (event  13  in  Table  3)  are  less  than  44  km 
in  depth.  The  one  exception  has  a  focal  depth  of  88  km.  For  some  events,  records 
from  several  stations  were  available,  whereas  for  some  events  only  a  single 
record  was  us'ed  because  of  the  poor  signal/noise  ratios  at  other  stations.  Depths 
identified  with  an  exclamation  mark  were  determined  by  Harvard  University  and 
those  identified  with  an  asterisk  were  obtained  in  this  study  by  matching  the 


9 


Table  2.  Earthquakes  studied  using  the  two-station  method 


Event  Date 

Origin 

(h:m:s) 

Lat 

(°N) 

Lon 

(°E) 

Depth 

(km) 

Recording 

stations 

Az-dev  1 
(°) 

Az-dev  2 
(°) 

Region 

1 

921006 

08:57:17.7 

38.42 

56.52 

10 

GNI-ANTO 

3.6 

1.9 

1 

2 

941129 

14:30:28.4 

38.71 

20.48 

21 

ANTO-GNI 

2.0 

1.0 

1 

3 

940329 

07:56:53.9 

29.10 

51.26 

33 

GNI-KIV 

7.6 

5.7 

1 

4 

940920 

05:51:46.0 

32.50 

48.77 

33 

GNI-KIV 

4.1 

2.7 

1 

5 

950422 

00:21:48.6 

30.89 

49.91 

25 

GNI-KIV 

5.8 

4.1 

1 

6 

940524 

02:05:36.2 

38.66 

26.54 

17 

GNI-ABKT 

7.1 

4.1 

2 

7 

940524 

02:18:34.9 

38.76 

26.60 

16 

GNI-ABKT 

7.8 

3.9 

2 

8 

940603 

22:40:21.2 

28.74 

70.07 

33 

ABKT.KIV 

6.9 

3.6 

2 

9 

931026 

10:03:57.5 

36.74 

28.05 

33 

ANTO-KIV 

6.3 

2.3 

2 

10 

950103 

22:51:45.0 

34.94 

23.55 

33 

ANTO-KIV 

1.3 

0.7 

2 

11 

950203 

22:29:13.3 

34.42 

25.03 

49 

ANTO-KIV 

6.7 

3.3 

2 

12 

950330 

18:17:15.1 

34.40 

24.79 

7 

ANTO-KIV 

5.8 

2.9 

2 

13 

960422 

14:42:32.4 

39.17 

47.37 

29 

AFIF-RANI 

7.2 

6.2 

3 

14 

951123 

18:07:17.3 

29.33 

34.75 

10 

AFIF-HALM 

4.7 

4.0 

3 

15 

951124 

16:43:45.5 

28.94 

34.71 

10 

AFIF-HALM 

6.8 

5.8 

3 

16 

951125 

11:41:35.3 

29.12 

34.84 

10 

AFIF-HALM 

5.5 

4.7 

3 

17 

951207 

17:48:16.3 

14.56 

55.77 

0 

HALM-AFIF 

7.6 

6.9 

3 

RAYN-BGIO 

3.0 

1.6 

3 

18 

960328 

07:28:28.1 

11.92 

57.81 

10 

HALM-AFIF 

4.9 

4.5 

3 

19 

960402 

07:59:26.2 

37.83 

27.00 

33 

AFIF-HALM 

6.6 

6.1 

3 

BGIO-HALM 

4.6 

1.8 

3 

20 

960426 

07:01:27.5 

36.37 

28.04 

75 

AFIF-HALM 

4,8 

4.5 

3 

21 

960225 

16:14:10.9 

35.65 

57.07 

33 

HALM-RANI 

5.2 

4.7 

3 

22 

960225 

17:42:04.5 

35.65 

57.05 

33 

HALM-RANI 

5.3 

4,7 

3 

23 

960331 

16:02:09.3 

31.83 

49.03 

0 

HALM-SODA 

1.8 

1.2 

3 

24 

960103 

08:42:26.3 

39.03 

48.72 

62 

RANI-SODA 

7.6 

6.6 

3 

25 

960128 

08:43:16.2 

34.26 

46.45 

33 

RANI-SODA 

6.1 

5,0 

3 

26 

940610 

03:00:44.5 

38.56 

70.63 

33 

ABKT-BGIO 

3.7 

1.3 

1,2,3 

27 

951026 

04:23:27.5 

39.18 

72.07 

50 

ABKT-BGIO 

0.4 

0.2 

1,2,3 

spectral  holes  of  the  theoretical  spectra  with  those  of  the  observed  data.  The  dip, 
slip  and  strike  angles  were  obtained  by  Harvard  University.  The  number  in  the 
rightmost  column  of  Table  3  indicates  the  region  that  the  smface  waves  traverse. 


5.  Interstation  Group  and  Phase  Velocities 

Before  using  the  two-station  method,  we  removed  the  mean  ground  motion 
from  the  records,  corrected  for  geometrical  spreading  and  instriunent  response, 
and  applied  bandpass  filtering.  Figure  4  is  an  example  showing  the  original  trace 
(top)  recorded  at  station  SODA  for  an  earthquake  that  occurred  on  28  January 
1996,  and  the  windowed  trace  of  ground  displacement  (bottom)  after  correcting 


10 


Table  3.  Earthquakes  studied  using  the  single-station  method 


Event 

Date 

Origin 

(h:m:s) 

Lat 

(°N) 

Lon 

(OE) 

Depth 

(km) 

Dip 

(°) 

Slip 

(°) 

Strike 

Recording 

stations 

Region 

1 

940920 

05:51:46.0 

32.50 

48.77 

15* 

25 

87 

103 

ABKT,  GNI 
BGIO 

1 

3 

2 

950422 

00:21:48.6 

30.89 

49.91 

8* 

32 

72 

117 

GNI 

1 

3 

940329 

07:56:53.9 

29.10 

51.26 

6* 

40 

104 

334 

ABKXKIV 

1 

4 

941120 

14:31:02.2 

35.34 

39.56 

29 

57 

176 

141 

KIV 

1 

5 

910724 

09:45:41.8 

36.52 

44.07 

10* 

44 

121 

335 

KIV 

1 

6 

940620 

09:09:02.9 

28.97 

52.61 

14* 

67 

-5 

251 

ANTO 

1 

7 

941214 

20:43:53.7 

35.10 

58.63 

5* 

32 

144 

319 

ANTO 

1 

8 

950223 

21:03:01.3 

35.05 

32.28 

15! 

21 

140 

239 

ANTO 

1 

9 

940701 

10:12:41.2 

40.23 

53.38 

44! 

30 

71 

251 

GNI,KIV 

2 

10 

940701 

19:50:04.3 

40.22 

53.39 

44 

34 

108 

290 

GNI, KIV 

2 

11 

910712 

10:42:21.2 

45.36 

21.06 

11 

78 

-1 

279 

KIV 

2 

12 

911202 

08:49:40.2 

45.50 

21.12 

9 

72 

-5 

103 

KIV 

2 

13 

900531 

00:17:47.8 

45.81 

26.77 

88 

26 

54 

90 

KTV 

2 

14 

960221 

04:59:51.2 

28.80 

34.78 

10 

30 

-104 

132 

RANI, SODA 

3 

15 

951123 

18:07:17.2 

29.33 

34.75 

26* 

77 

7 

199 

APIP,HALM 

3 

16 

960226 

08:08:19.2 

28.32 

57.09 

33 

7 

125 

315 

RANI, SODA 

3 

17 

960314 

21:47:57.9 

14.74 

55.74 

18* 

20 

-118 

72 

AFIF,RANI, 

HALM,SODA 

3 

18 

951205 

18:49:30.4 

39.44 

40.15 

29! 

70 

160 

136 

BGIO 

1,3 

for  geometrical  spreading  and  instrument  response,  and  bandpass  filtering 
between  comer  frequencies  of  0.01  and  1  Hz.  The  epicentral  distance  and  the 
starting  times  of  the  traces  from  the  origin  are  indicated  in  the  boxes.  The  dashed 
lines  above  the  top  trace  represent  the  window  with  a  10%  cosine  taper.  In  this 
example,  the  window  length  is  400  sec,  and  the  starting  and  ending  times  corre¬ 
spond  to  group  velocities  of  about  2.2  to  4.5  km/sec,  respectively. 

The  Multiple  Filter  Technique  (MFT)  was  then  utilized  to  extract  group 
velocities.  Figrure  5  shows  contoured  spectral  amplitudes  of  fundamental-  and 
higher-mode  Rayleigh  waves  obtained  for  an  earthquake  that  occurred  on  22 
April  1995  and  recorded  at  station  GNI.  The  windowed  ground  displacement, 
after  correction  for  the  geometrical  spreading  and  instrument  response,  and 
bandpass  filtering,  is  displayed  at  the  rightmost  part  of  the  figure.  The  same 
trace  with  a  linear  group  velocity  scale  is  plotted  adjacent  to  the  dispersion  dia¬ 
gram.  This  trace  can  be  used  to  correlate  group  velocity  contour  values  to  the 
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trace  itself.  The  long,  smooth  curve  composed  of  squares  identifies  the  fundamen¬ 
tal-mode  group  velocity  dispersion.  The  shorter  curve  with  a  group  velocity  of 
about  3.1  km/s  at  5  sec  is  a  superposition  of  higher  modes.  The  square,  circle,  tri¬ 
angle,  and  plus  signs  correspond  to  the  fom-  largest  amplitudes  in  the  time 
domain  when  the  MFT  is  applied. 

A  phase-matched  filter  technique  as  proposed  by  Herrin  and  Goforth  (1977), 
and  coded  by  Herrmann  (1987)  was  used  to  extract  the  phase-velocity  informa¬ 
tion.  The  interstation  group  and  phase  velocities  can  then  be  obtained  for  the 
three  regions  as  shown  in  Figures  6,  7,  and  8.  In  each  of  the  three  figures,  group 
velocities  appear  below  the  phase  velocities. 

Different  features  can  be  observed  for  the  dispersion  curves  of  the  three 
regions.  For  Region  1  (Figure  6),  the  group  velocities  obtained  for  different  events 
are  similar  to  one  another  at  periods  below  20  sec  and  display  greater  scatter  at 
longer  periods.  The  velocities  increase  rapidly  with  increasing  period  at  shorter 
periods,  and  reach  a  value  of  about  2.89  km/s  at  10  sec.  The  velocity  values 
remain  almost  unchanged  in  the  period  range  10-20  sec.  Above  20  sec,  the  veloci¬ 
ties  increase  again,  smoothly  and  monotonically,  reaching  their  maximum  values 
at  the  maximum  observable  period.  The  phase  velocity  curves  for  this  region, 
however,  increase  uniformly  throughout  the  entire  period  range. 

Most  of  the  curves  for  Region  2  (Figure  7)  are  similar  to  one  another,  but 
velocities  in  a  few  cases  deviate  substantially  from  mean  values.  This  deviation 
may  be  due  to  systematic  errors  caused  by  lateral  variations  of  properties  in  this 
region.  It  is  comprised  of  both  oceanic  (Black  and  Caspian  Seas)  and  continental 
portions.  Both  phase  and  group  velocities,  at  shorter  periods,  are  lower  than 
those  for  Region  1.  This  implies  that  the  shear  velocities  at  shallow  depths  in 
Region  2  are  lower  than  those  in  Region  1. 

Fifteen  interstation  group  and  phase  velocities  are  obtained  for  Region  3 
(Figure  8).  The  largest  interstation  distance  is  1370  km  between  stations  BGIO 
and  RAYN.  The  distances  for  the  station  pairs  AFIF-RANI,  AFIF-HALM,  HALM- 
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RANI,  and  HALM-SODA  are,  however,  in  the  range  178-543  km.  Short  intersta¬ 
tion  distances  can  cause  problems  when  using  the  two-station  method.  Any 
departures  from  a  great-circle  path  or  measurement  errors  can  result  in  severe 
bias  in  the  velocity  and  attenuation  determinations.  To  control  the  errors,  we 
tried  to  collect  as  many  events  for  each  station  pair  as  possible.  For  example, 
seven  events  were  used  for  the  station  pair  AFIF-HALM.  If  observed  values  are 
similar,  the  standard  deviations  of  the  mean  values  will  be  reduced.  In  addition, 
results  for  the  station  pairs  with  shorter  interstation  distances  will  be  compared 
with  those  with  longer  interstation  distances,  and  also  be  compared  with  the 
average  values  to  examine  the  usability  of  the  data.  The  curves  for  the  station 
pairs  RAYN-BGIO  and  BGIO-HALM  lie  in  the  middle  of  the  group,  suggesting 
that  the  selected  cmwes  for  the  shorter  interstation  distances  are  realistic.  The 
models  obtained  after  inversion  will  be  compared  with  those  reported  in  the  pre¬ 
vious  studies  to  check  the  validity  of  the  interstation  results.  The  mean  value  of 
group  velocities  at  10  sec  for  Region  3  is  higher  than  that  foxmd  for  either  Region 
1  or  2,  being  about  2.98  km/s.  This  suggests  that  the  near-surface  shear  velocity 
is  higher  here  than  in  the  other  two  regions. 

The  interstation  group  and  phase  velocities  for  the  station  pair  ABKT-BGIO 
are  plotted  in  Figure  9.  The  paths  connecting  the  two  stations  cross  the  southern 
edge  of  the  Caspian  Sea,  the  Turkish  and  Iranian  Plateaus,  and  the  northern 
Arabian  Peninsula  as  shown  in  Figure  2.  According  to  our  regionalization,  the 
paths  cross  all  of  the  three  regions.  Since  our  goal  is  to  obtain  velocity  and  atten¬ 
uation  structure  for  individual  regions,  we  did  not  use  the  dispersion  data  for  this 
station  pair  in  our  inversion.  If  our  regionalization  is  appropriate,  i.e.  the  three 
regions  are  geologically  different  with  different  velocity  and  attenuation  struc¬ 
ture,  the  group  and  phase  velocities,  and  the  attenuation  values  for  this  station 
pair  should  give  an  average  for  the  three  regions.  The  average  group  velocities  at 
10  sec  are  about  2.89,  2.35,  and  2.98  km/s  for  Regions  1,  2,  and  3,  respectively. 
The  average  group  velocity  at  10  sec  for  these  two  events  is  about  2.70  km/s. 
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which  is  in  the  middle  of  the  three  values.  The  velocities  increase  with  period  and 
have  a  value  of  about  3.38  km/s  at  40  sec,  which  is  also  in  the  middle  of  the  aver¬ 
age  values  of  3.33,  3.30,  and  3.55  km/s  for  Regions  1,  2,  and  3  at  that  period, 
respectively.  The  group  velocities  for  these  two  events  are  higher  than  those  for 
Region  2,  lower  than  those  for  Region  3,  and  similar  to  those  for  Region  1  for  the 
periods  greater  than  20  sec.  The  average  phase  velocity  for  these  two  events  at  10 
sec  is  about  3.12  km/s.  This  value  is  also  in  the  middle  of  the  corresponding  val¬ 
ues,  3.13,  2.82,  and  3.28  for  the  three  regions.  From  the  above  comparisons,  we 
conclude  that  our  division  of  the  Middle  East  into  three  regions  is  appropriate. 

6.  Shear  Velocity  Models  Obtained  Using  the  Two-station  Method 

Since  the  dispersion  curves  in  each  group  exhibit  similar  features,  the  aver¬ 
age  values  and  the  standard  errors  were  calculated  for  each  of  the  three  regions. 
It  should  be  noted,  however,  that  for  some  periods  only  single  observations  are 
available  and  the  standard  deviations  cannot  be  calculated.  In  these  cases,  we 
compared  the  measurement  errors  at  these  periods  with  the  estimated  standard 
deviations  at  other  periods  where  multiple  observations  are  available,  choosing 
the  largest  values  as  our  conservative  error  estimates.  We  chose  0.3  km/s  as  the 
error  estimates  at  those  points. 

To  save  computing  time,  the  average  dispersion  data  were  simplified  to 
make  a  roughly  even  distribution  over  a  linear  period  range.  The  simplified 
group  and  phase  velocities  were  then  inverted  using  a  dififerential  inverse  tech¬ 
nique  to  obtain  the  shear- velocity  models  for  the  three  regions.  In  the  inversion 
process,  the  thickness  of  each  layer  was  fixed.  The  initial  values  of  compressional 
and  shear  velocities,  and  densities  are  6.0  km/s,  3.5  km/s,  and  3.0  g/cm^  for  all  of 
the  layers.  These  values  were  changed  later  to  see  the  effect  of  initial  values  on 
inversion  results.  Figure  10  is  an  example  showing  the  shear  velocity  models  and 
the  resolving  kernels  obtained  by  inverting  the  fimdamental-mode  Rayleigh  data 
for  Region  1.  The  left  panel  of  the  figure  shows  the  velocity  model,  whereas  the 
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Figure 


right  panel  gives  the  normalized  resolving  kernels  at  six  depths.  The  dashed  line 
in  each  figure  represents  the  initial  value  of  the  shear  velocities,  3.5  km/sec.  The 
lateral  bars  sire  the  error  estimates  of  the  corresponding  velocity  values,  plotted 
at  the  midpoint  of  the  corresponding  layers.  Figure  11  shows  the  observed  group 
and  phase  velocities  as  well  as  the  theoretical  dispersion  curves  predicted  by  the 
model  obtained  for  Region  1.  The  circles  denote  the  group  velocities,  whereas  the 
squares  represent  the  phase  velocities.  The  vertical  bars  are  the  standard  errors 
of  the  data.  The  curves  for  the  resolving  kernels  on  the  right  panel  of  Figure  10 
show  that  the  shallow  structure  is  better  resolved  tham  the  deeper  structure.  The 
shapes  of  the  velodty  resolving  kernels  are  similar  for  all  of  the  three  regions, 
and  are  therefore  not  plotted  for  the  other  two  regions.  The  velocity  models 
obtained  in  this  study,  from  analysis  of  the  resolving  kernels,  are  restricted  to  the 
crust  and  the  uppermost  part  of  mantle  within  a  depth  of  about  120  km.  The 
match  of  the  theoretical  velocity  dispersion  curves  to  the  observed  data  (Figure 
11)  is  good  over  the  entire  period  range.  Similar  results  were  obtained  for  the 
other  two  regions,  and  are  not  plotted.  For  comparison,  we  plot  the  shear- velocity 
models  for  the  three  regions  together  in  Figure  12.  The  dotted,  dashed,  and  solid 
lines  are  the  models  for  Regions  1,  2,  and  3,  respectively.  The  uncertainties  of  the 
velocities  in  each  layer  are  also  plotted  as  horizontal  bars,  which  overlap  for  the 
deeper  layers.  The  velocities  for  the  three  regions  differ  greatly  in  the  upper  10 
km.  Region  3  has  the  highest  velocity  values,  whereas  Region  2  has  the  lowest 
velocities.  The  average  velocities  for  the  upper  10  km  are  about  3.21,  2.85,  and 
3.39  km/s  for  Regions  1,  2,  and  3,  respectively.  Below  10  km  in  the  crust,  the 
velocities  in  Regions  2  and  3  increase  gradationally  without  any  obvious  gradient 
change,  whereas  the  velocity  gradient  in  Region  1  changes  significantly  at  depths 
between  about  12  and  30  km.  In  the  uppermost  mantle,  the  velocity  gradients 
decrease  for  all  of  the  three  regions,  and  velocity  values  in  the  uppermost  mantle 
do  not  differ  much  fi:om  region  to  region,  especially  for  Regions  1  and  3.  At  depths 
near  120  km  the  velocities  for  the  three  regions  approach  a  high  value  in  the 
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Figure  12 


range  4.64-4.78  km/s.  No  low-velocity  zones  are  resolved  for  any  of  the  three 
regions  in  the  present  study.  Since  Region  2  includes  both  continental  and 
oceanic  portions,  the  model  for  this  region  may  not  represent  either  of  those  por¬ 
tions  very  well. 

Figure  13  compares  our  velocity  models  with  those  of  other  studies  in  the 
Arabian  Peninsula.  The  bold  line  is  the  model  obtained  in  this  study,  whereas  the 
thin  lines  and  dashed  or  dotted  lines  represent  the  results  from  other  studies. 
The  models  labeled  by  "Crust  H"  and  "Crust  L"  are  cited  from  the  study  of 
Knopoff  and  Fouda  (1975).  Their  study  showed  a  variety  of  solutions  which  can 
explain  the  data.  We  selected  two  typical  models,  in  the  depth  range  0-180  km, 
from  the  Crust  H  and  Crust  L  models  to  compare  with  our  model. 

Mokhtar  and  Al-Saeed  (1994),  using  data  recorded  by  station  RYD  in  the 
center  of  Arabian  Peninsula,  studied  three  major  propagation  paths.  Path  I  con¬ 
nects  the  station  and  the  events  in  the  Red  Sea,  Path  II  is  for  the  events  in  Gulf 
of  Aden,  and  Path  III  is  for  the  events  in  the  Zagros  folded  zone.  The  results  are 
plotted  in  Figure  13  using  the  lines  labeled  by  Path  I,  Path  II,  and  Path  III, 
respectively.  Because  the  period  range  of  the  dispersion  data  is  from  4  to  42  sec, 
the  maximum  depth  for  the  models  is  60  km. 

Ghalib  (1992)  used  data  recorded  at  several  stations  surrounding  the  Ara¬ 
bian  Peninsula  to  invert  for  the  velocity  structure  using  a  single-station  method. 
Since  his  maximum  period  is  about  70  sec,  the  maximum  depth  he  used  for  the 
inversion  is  100  km.  Many  paths  were  included  in  his  study.  For  purposes  of  com¬ 
parison,  we  select  only  two  typical  paths  crossing  the  entire  Arabian  Peninsula, 
one  trending  north-south,  another  trending  east-west.  The  models  for  the  two 
paths  are  plotted  in  Figure  13  as  the  lines  labeled  by  T190  and  S212f. 

Figure  13  shows  that  the  shear-velocity  model  obtained  in  this  study  is  in 
good  agreement  with  models  obtained  in  other  studies,  especially  for  the  depth 
range  20-40  km.  For  the  upper  20  km,  our  velocities  agree  with  the  average  of 
other  studies  because  we  used  paths  at  several  azimuths,  and  averaged  the 
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dispersion  values  to  obtain  an  average  velocity  model  for  the  entire  peninsula. 
For  the  uppermost  mantle,  at  the  depths  40-100  km,  there  is  a  wide  range  of 
velocities  (4.4-4. 7  km/s).  The  model  obtained  in  this  study  (4.4  km/s)  lies  in  the 
lower  portion  of  that  region. 

The  velocity  model  of  Region  1  is  compared  with  those  of  Mindevalh  and 
Mitchell  (1989)  as  shown  in  Figure  14.  In  their  study,  the  region  is  divided  into 
two  subregions;  eastern  and  western  Turkey.  The  figure  shows  that  the  velocity 
model  of  this  study  is  in  good  agreement  with  their  models.  In  the  lower  crust 
and  the  uppermost  mantle  (20-80  km)  the  three  models  are  nearly  identical.  For 
the  upper  crust,  the  model  of  this  study  agrees  well  with  the  model  for  eastern 
Turkey,  and  is  slightly  higher  than  that  for  western  Turkey.  This  is  not  surprising 
because  the  path  connecting  stations  ANTO  and  GNI  is  located  within  eastern 
Turkey.  In  addition  to  using  data  for  the  station  pair  ANTO-GNI,  we  have  used 
data  from  the  station  pair  GNI-KIV  to  determine  an  average  velocity  model  for 
this  region.  The  model,  therefore,  appears  to  be  valid  for  the  motmtain  area 
between  the  Black  Sea  and  Caspian  Sea  as  well  as  for  eastern  Turkey. 

In  the  inversion  process,  we  used  an  initial  velocity  value  of  3.5  km/s  for  all 
layers.  The  velocity  models  obtained  can  be  affected  by  the  initial  model  values, 
especially  at  depths  where  the  resolving  kernels  are  broad.  To  investigate  the 
effect  of  initial  values  on  the  inversion  results,  we  performed  the  inversion  two 
more  times,  using  initial  values  of  4.0  and  4.5  km/s  separately  for  all  layers.  The 
predicted  group  and  phase  velocities,  and  the  resolving  kernels  are  almost  identi¬ 
cal  to  those  previously  obtained.  They  are,  therefore,  not  plotted  here.  The  veloc¬ 
ity  models  obtained  using  different  initial  values  for  the  three  regions  are  pre¬ 
sented  in  Figure  15.  The  bold,  thin,  and  dashed  lines  in  each  set  of  curves  corre¬ 
spond  to  the  initial  values  of  3.5,  4.0,  and  4.5  km/s,  respectively. 

Figxire  15  indicates  that  the  velocities  at  crustal  depths  are  not  changed  by 
the  variation  of  the  initial  values.  At  mantle  depths  the  changes  in  velocity  are 
very  limited,  far  less  than  the  error  estimates.  None  of  the  curves  approaches  the 
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initial  values  at  great  depths.  We  conclude  that  the  models  for  the  three  regions 
are  not  affected  by  the  initial  values  throughout  the  depth  range.  Conservatively 
speaking,  the  velocities  can  be  resolved  to  at  least  120  km,  a  depth  similar  to  that 
suggested  by  the  resolving  kernels. 

7.  Shear-wave  Q  Models  Obtained  Using  the  Two-station  Method 

Figure  16  shows  the  interstation  attenuation  coefficients  (y)  for  Region  3. 
As  with  velocities,  we  calculated  the  average  value  of  y  for  each  of  the  three 
regions.  For  the  periods  where  only  a  single  measurement  is  available,  a  value  of 
0.5x10"^  km  ^  was  chosen  as  a  conservative  error  estimate.  The  average  attenua¬ 
tion  data  were  then  simplified  to  make  a  roughly  even  distribution  over  a  linear 
period  range.  The  differential  inverse  procedure  was  used  to  obtain  the  mod¬ 
els  for  the  three  regions,  taking  0.0  as  an  initial  value  of  for  all  layers.  Fig¬ 
ure  17  is  an  example  showing  the  models  and  the  resolving  kernels  for  Region 
3.  The  final  model  with  standard  deviations  is  plotted  along  with  the  starting 
model  in  the  left  panel.  The  right  panel  gives  the  normalized  resolving  kernels  at 
six  depths.  The  resolving  kernels  for  are  similar  for  all  of  the  three  regions, 
and  those  for  the  other  two  regions  are,  therefore,  not  plotted. 

Figure  17  shows  that  the  resolving  kernels  for  are  wider  at  great  depths 
than  those  for  the  sheeu:  velocity.  This  indicates  that  the  applicable  depth  of  our 
Qf,  models  are  shallower  than  that  of  shear  velocity.  We  restrict  the  models  of 
this  study  to  the  uppermost  80  km.  The  observed  attenuation  coefficients  and 
theoretical  dispersion  curves  predicted  by  the  model  for  Region  3  are  plotted  in 
Figure  18.  The  agreement  between  them  is  good  over  the  entire  period  range. 
Similar  results  were  obtained  for  the  other  two  regions,  and  are  not  plotted.  The 
models  obtained  for  the  three  regions  are  presented  in  Figure  19.  The  dotted, 
dashed,  and  solid  lines  are  the  models  for  Regions  1,  2,  and  3,  respectively.  The 
uncertainties  in  values  for  each  layer  are  also  plotted  as  horizontal  bars, 
which  overlap  for  deeper  layers.  The  models  show  different  features  from 
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region  to  region.  in  Region  3  is  highest  through  the  upper  60  km,  whereas 
in  Region  1  is  the  lowest  in  the  same  depth  range.  For  the  upper  10  km,  the  aver¬ 
age  Qf,  values  for  Regions  1,  2,  and  3  are  63,  71,  and  201,  respectively.  values 
for  Regions  1  and  3  decrease  with  depth  in  the  upper  crust,  whereas  for 
Region  2  increases.  The  model  for  Region  1  has  a  minimxim  at  a  depth  of  20 
km.  The  model  for  Region  3  exhibits  a  zone  of  constant  Q  in  the  depth  range 
20-80  km.  The  model  for  Region  2  fluctuates,  having  a  minimum  and  a  maxi¬ 
mum  at  depths  of  about  12  and  30  km,  respectively.  This  feature  may  be 
specious,  and  is  probably  due  to  systematic  errors  in  our  one-dimensional  inver¬ 
sion  in  a  region  where  varies  rapidly  in  three  dimensions.  As  with  the  discus¬ 
sion  for  velocities,  we  think  that  details  of  the  model  for  this  region  should  be 
ignored.  At  depths  greater  than  80  km,  the  values  for  all  three  regions 
increase  gradually,  approaching  the  initial  value. 

We  originally  used  a  initial  value  of  0.0  for  Qf  in  the  inversion.  We  tested 
the  effect  of  initial  values  of  on  the  resulting  models  by  inverting  using  ini¬ 
tial  values  of  0.005  and  0.01  to  (corresponding  to  values  of  200  and  100). 
The  resulting  models  for  the  three  regions  are  presented  in  Figure  20.  The 
models  are  more  easily  affected  by  the  initial  values  than  are  the  velocity  models, 
indicating  that  the  resolvable  depth  is  less  than  that  for  the  velocity.  At  crustal 
depths,  none  of  the  models  for  the  three  regions  is  greatly  affected  by  the  varia¬ 
tion  of  initial  values,  indicating  that  resolution  is  good  at  those  depths.  The 
depths  at  which  the  three  curves  corresponding  to  the  three  initial  values  begin 
to  deviate  are  about  80,  120,  and  90,  for  Regions  1,  2,  and  3,  respectively.  At 
greater  depths,  the  values  approach  the  initial  values,  indicating  poor  resolu¬ 
tion  at  those  depths.  To  be  conservative,  we  choose  80  km  as  the  applicable  depth 
of  the  Qf,  models  for  all  three  regions.  All  of  our  conclusions  concerning  the 
models  will  be  restricted  to  this  depth.  This  is  also  in  agreement  with  the  resrdt 
obtained  by  analysis  of  the  resolving  kernels.  Our  results  indicate  that  the 
resolvable  depth  for  Q^,  when  inverting  Rayleigh  wave  data,  is  less  than  that  for 
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Figure  20 


shear  velocity. 


8.  Shear-wave  Q  Models  Obtained  Using  the  Single-station  Method 

Determination  of  surface  wave  attenuation,  using  the  two-station  method, 
requires  that  the  azimuths  between  earthquakes  and  stations  are  similar  to  the 
azimuths  of  the  line  connecting  the  two  stations.  This  greatly  reduces  the  number 
of  usable  earthquakes.  In  addition,  if  the  two  stations  are  located  in  different  geo¬ 
logical  provinces,  only  the  average  attenuation  of  the  different  provinces  between 
the  two  stations  can  be  obtained.  For  these  reasons,  a  single-station  method  was 
developed  by  Cheng  and  Mitchell  (1981)  for  the  pxn-pose  of  improved  regional 
attenuation  studies.  The  importance  of  the  method  lies  in  its  ability  to  utilize  rel¬ 
atively  short  paths  between  a  single  source  and  a  single  receiver,  both  of  which 
are  preferably  located  within  a  single  geological  province. 

The  method  attempts  to  match  theoretical  and  observed  surface  wave  spec¬ 
tra.  This  is  done  by  a  trial-and-error  procedure  in  this  study.  The  method  utilizes 
both  fundamental-mode  and  higher-mode  spectra  since  use  of  the  amplitude 
ratios  for  the  two  spectra  obviates  the  need  to  know  the  seismic  moment  of  the 
earthquake  used.  To  extract  the  observed  spectnnn,  the  multiple-filter  technique 
was  applied  to  the  data  as  described  in  the  previous  section.  Computations  of  the¬ 
oretical  surface-wave  spectra  requires  knowledge  of  the  source  mechanism  of  the 
earthquake.  For  all  of  the  earthquakes  we  used  the  focal  mechanism  determined 
by  Harvard  University. 

Figure  21  is  an  example  showing  the  process  of  determination  for  each 
layer.  The  earthquake  in  this  example  occurred  on  March  29,  1994  (event  3  in 
Table  3)  as  recorded  at  two  stations  KIV  and  ABKT.  Station  names  are  shown  to 
the  left  of  the  corresponding  spectra.  In  the  tests,  three  layers  were  assumed  with 
layer  thicknesses  of  10,  30,  and  40  km,  overlying  a  halfspace.  The  values  in 
the  layers  are  denoted  by  Ql,  Q2,  Q3,  and  Q4.  In  the  first  test  (Figure  21a),  Q1 
varies  between  30  and  130,  with  a  step  of  20.  Q2,  Q3,  and  Q4  are  fixed  at  50,  90, 
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and  700,  respectively.  These  values  were  inferred  from  the  results  obtained  using 
the  two-station  method.  The  Vcdue  in  the  halfspace,  Q4,  is  somewhat  arbitrary 
since  we  already  know  that  cannot  be  resolved  at  those  depths.  As  can  be  seen 
by  comparing  the  spectra  for  both  stations,  the  only  acceptable  solution  for  Ql, 
among  the  values  shown,  is  70.  When  Ql  is  90,  the  fit  for  the  higher  modes  is 
good,  but  at  shorter  periods  theoretical  values  of  the  fundamental  mode  for  sta¬ 
tion  KIV  are  too  high. 

Figure  21b  shows  comparisons  of  the  spectra  when  Q2  varies.  In  this  test, 
Ql  takes  the  value  of  70  as  obtained  in  the  previous  test.  Q3  and  Q4  were  fixed  at 
90  and  700,  respectively.  Q2  varies  between  10  and  110  with  a  step  of  20.  Both 
fimdamental-  and  higher-mode  spectra  are  sensitive  to  the  variation  of  Q2.  When 
Q2  is  less  than  50,  the  theoretical  cruwes  are  lower  than  observed  values  for  both 
fundamental  amd  higher  modes  at  the  two  stations.  When  Q2  is  greater  than  50, 
the  opposite  occurs.  The  only  good  match  occurs  when  Q2  is  50. 

Figure  21c  shows  the  comparison  when  Q3  varies.  In  this  case  Ql  and  Q2 
are  fixed  at  70  and  50,  respectively,  as  determined  before,  and  Q4  is  fixed  at  700. 
Q3  varies  between  30  and  180  with  a  step  of  30.  The  theoretical  curves  for  the 
fundamental  mode  are  far  less  sensitive  to  the  variation  of  Q3  than  are  the 
higher  modes.  A  good  match  can  be  detected  when  Q3  takes  the  values  between 
60  and  120.  A  middle  value  of  90  is  therefore  chosen  as  the  solution  for  Q3. 

In  Figure  21d  we  assign  70,  50  and  90  to  Ql,  Q2,  and  Q3,  allowing  Q4  to 
vary  from  400  to  2000.  No  spectral  changes  can  be  detected  when  Q4  varies.  This 
indicates  that  the  spectra  are  not  sensitive  to  changes  in  at  great  depths.  No 
sensitivity  to  at  depths  greater  than  80  km  was  found  for  any  case  tested. 
Therefore,  the  models  obtained  using  the  single-station  method  are  restricted 
to  the  upper  80  km  in  this  study,  the  same  depth  range  determined  in  data  pro¬ 
cessing  using  the  two-station  method. 

The  Qft  models  for  the  three  regions  obtained  using  the  single-station 
method  are  listed  in  Table  4,  and  superposed  on  those  obtained  using  the  two- 
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station  method  (bold  lines)  as  shown  in  Figure  22.  Comparisons  of  theoretical 
spectra  predicted  by  the  corresponding  models  with  the  observed  data  for  all  the 
earthquakes  studied  in  Regions  1,  2,  and  3  are  presented  in  Figures  23,  24,  and 
25,  respectively.  The  events,  the  recording  stations,  and  models  obtained  are 
indicated  in  the  lower  right  comers  of  the  panels. 


Table  4.  models  obtained  using  the  single-station  method 


Layer 

Depth 

Q. 

Qft 

Q. 

Q. 

Number 

(km) 

(model  1) 

(model  2) 

(model  3) 

(model  4) 

Region  1 

1 

0-10 

60 

70 

60 

50 

2 

10-40 

50 

50 

60 

50 

3 

40-80 

90 

90 

90 

90 

4 

>80 

700 

Region  2 

700 

700 

700 

1 

0-10 

60 

2 

10-20 

100 

3 

20-40 

70 

4 

40-80 

120 

5 

>80 

700 

Region  3 

1 

0-20 

180 

140 

40  (for  0  —  2  km) 
180  (for  2  -  20  km) 

2 

20-80 

140 

140 

140 

3 

>80 

700 

700 

700 

For  events  1  and  3  in  Region  1  (Figure  23),  two  seismograms  are  available. 
In  general,  the  match  is  good  for  both  fundamental-mode  and  higher-mode  spec¬ 
tra  for  all  of  the  events  used.  For  this  region,  four  models  were  obtained,  indicat¬ 
ing  that  the  models  are  path  dependent.  The  fom  models  differ  only  in  the 
cmst.  varies  between  50  and  70  for  the  upper  10  km,  and  between  50  and  60 
for  the  xmderl3dng  layer  (10-40  km).  Although  the  four  models  differ,  they  all  lie 
close  to  the  model  obtained  by  the  two-station  method  (bold  line).  This  indicates 
that  the  model  obtained  using  the  two-station  method  for  this  region  is  applicable 
and  can  be  thought  of  as  an  average  model  for  the  region. 

One  of  the  advantages  of  using  the  single-station  method  is  that  the  appli¬ 
cable  range  for  the  models  obtained  using  the  two-station  method  can  be 
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Figure  22 
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Figure  24 
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broadened.  We  previously  used  the  station  pairs  ANTO-GNI  and  GNI-KIV  in 
Region  1.  The  results  for  velocity  and  are  only  valid  for  the  Turkish  Plateau 
or,  strictly  speaking,  the  region  along  the  interstation  paths.  Figure  3  shows  that 
some  paths  lie  in  the  Iranian  Plateau  and  some  paths  go  through  both  the  Turk¬ 
ish  and  Iranian  Plateaus.  This  leads  us  to  conclude  that  the  velocity  and 
models  previously  obtained  for  the  Turkish  Plateau  also  apply  to  the  Iranian 
Plateau. 

Comparisons  of  the  theoretical  and  observed  spectra  for  Region  2  are  plot¬ 
ted  in  Figure  24.  Five  events  recorded  at  stations  GNI  and  KIV  were  used.  Three 
paths  cross  the  northern  portion  of  the  Black  Sea,  whereas  4  paths  cross  the 
Caspian  Sea.  Events  9  and  10  are  recorded  at  both  stations.  All  paths  include 
both  oceanic  and  continental  segments. 

The  method  used  for  Region  1  is  not  applicable  to  Region  2  since  no  higher¬ 
mode  spectra  were  observable  for  any  of  the  5  events.  In  this  case,  we  first  simpli¬ 
fied  the  previously  obtained  model  for  this  region,  and  computed  only  the  funda¬ 
mental-mode  spectrum.  Strictly  speaking,  the  model  listed  in  Table  4  for  this 
region,  therefore,  is  not  the  one  obtained  using  the  single-station  method. 

Higher-mode  energy  is  probably  absent  because  it  cannot  propagate 
through  the  oceanic  cnist  beneath  the  Caspian  and  Black  Seas.  The  model 
obtained  for  this  region  is  less  imique  than  models  obtained  using  multimode 
spectra.  Fluctuations  in  the  model  may  not  represent  a  real  feature  of  the  crust, 
but  may  be  due  to  systematic  errors  in  the  attenuation  data  caused  by  the  lateral 
variations  of  and  velocity  in  this  region. 

Five  events  recorded  at  5  stations  were  used  for  Region  3  (Figure  25). 
Events  14  and  15  are  located  at  the  southern  end  of  the  Dead  Sea  fault.  Paths 
connecting  these  two  events  and  stations  cross  the  central  and  western  portions 
of  the  peninsula.  Event  17  lies  in  the  Gulf  of  Aden  and  the  path  to  the  station  tra¬ 
verses  small  portion  of  the  Gulf.  Most  of  the  paths  are  in  the  central  and  south¬ 
ern  part  of  the  peninsula.  Event  16  is  in  the  Zagros  folded  zone.  The  path  from 
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that  event  crosses  the  shield,  the  shelf,  and  the  basins. 

Panels  1  and  2  are  for  event  14  recorded  at  stations  RANI  and  SODA, 
respectively.  For  this  event,  model  1  (180  for  the  upper  20  km  and  140  for  depth 
range  20-80  km)  was  obtained  for  both  of  the  two  paths.  Panels  3  and  4  are  for 
event  15,  which  is  located  very  close  to  event  14.  For  this  event,  the  records  of 
two  stations  (AFIF  and  HALM)  are  available  for  modehng.  Since  the  paths 
deviate  only  slightly  from  that  connecting  event  14  and  RANI,  we  tried  the  same 
model  as  for  event  14.  As  shown,  the  theoretical  spectrum  for  the  fundamental 
mode  is  in  good  agreement  with  the  observed  data  at  each  station.  For  higher 
modes,  however,  the  observed  spectra  are  higher  than  expected  by  nearly  an 
order  of  magnitude.  Even  using  an  extremely  high  value  of  1000  for  all  of  the 
layers,  the  synthetics  for  the  higher  modes  are  still  lower  than  the  observed  spec¬ 
tra.  This  suggests  that  the  large  amplitudes  of  Lg  in  these  two  records  ceumot  be 
explained  by  high  in  the  crust.  Focusing  or  site  effects  that  affect  high-fre¬ 
quency  waves  coming  from  northwest  of  the  peninsula  might  be  responsible  for 
the  unusual  high  values  of  the  spectra. 

The  paths  from  event  16  to  the  recording  stations  SODA  and  RANI  cross 
basins  where  a  thick  cover  of  low-Q  sediment  exists.  for  those  paths  (model  2 
in  Table  4)  is  140  for  the  upper  20  km  of  the  crust.  This  value  is  22%  lower  than 
for  model  1.  Panels  7  and  8  in  Figure  25  show  an  optional  model  (model  3)  for 
the  same  paths  as  panels  5  and  6.  In  this  model,  in  the  uppermost  2  km  is  40 
and  in  the  next  18  km  is  180  as  in  model  1.  A  model  with  a  2-km  thick  layer  with 
a  of  40  can  explain  the  observed  spectra  in  the  eastern  portion  of  the  Arabian 
Peninsula.  Alternative  models  can  also  explain  the  spectra.  One  possibility  is  a 
3-km  thick  layer  in  which  is  50  and  another  is  a  4-km  thick  layer  in  which 
is  60.  The  predicted  spectra  are  nearly  identical  to  those  in  panels  7  and  8,  and 
are,  therefore,  not  plotted. 

Panels  9  through  12  are  for  event  17  and  the  recording  stations  AFIF, 
RANI,  HALM,  and  SODA,  respectively  For  those  paths,  model  2  provides  the 
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best  fit  to  the  observed  spectra.  Since  about  20-25%  of  those  paths  cross  water,  we 
might  expect  that  the  best  model  for  the  southern  portion  of  the  peninsula  would 
be  model  1  rather  than  model  2.  The  observed  spectra,  however,  exhibit  small 
spectral  holes  at  about  12  sec  that  are  not  predicted  by  the  given  focal  mecha¬ 
nism.  This  may  mean  that  the  reported  focal  mechanism  is  incorrect  and  our 
values  at  shallow  depths  may  be  biased  too  low  for  this  reason. 

In  general,  for  the  western,  central,  and  southern  portions  of  the  Arabian 
Peninsula,  the  crust  and  uppermost  mantle  can  be  modeled  as  a  20-km  thick 
layer  with  of  180  and  a  60-km  thick  layer  with  of  140.  For  the  eastern  por¬ 
tion  of  the  peninsula,  the  low  Q  sediments  must  be  taken  into  account.  for  the 
sediments  can  be  as  low  as  40-60. 

For  Regions  1  and  3,  the  models  obtained  using  the  two-  and  single-sta¬ 
tion  methods  show  good  agreement.  For  Region  2,  however,  models  are  not 
verified  by  the  single-station  method.  Since  the  higher-mode  energy  was  not 
detected  in  the  observation,  we  could  not  vary  the  value  of  any  layer,  but  used 
the  simplified  model  to  compute  the  theoretical  fundamental-mode  spectrum  and 
compared  it  to  the  observed  spectrum.  The  results  are  probably  less  unique  than 
models  obtained  using  multimode  spectra.  For  both  methods  the  paths  cross  both 
continental  and  oceanic  areas,  that  are  charactered  by  different  seismic  velocity 
and  distributions.  We  attempted  to  model  this  laterally  complex  region  using 
one-dimension  inversion.  We  think  that  the  fluctuations  in  the  model  may  be 
caused  by  the  lateral  complexities  of  and  velocity,  and  may  not  represent  the 
real  crustal  properties  of  this  region. 

Our  results  indicate  that  is  higher  in  the  Arabian  Peninsula  than  in  the 
Turkish  and  Iranian  Plateaus.  This  is  consistent  with  the  results  obtained  from 
Lg  coda  Q  measurements  (Cong  and  Mitchell,  1998).  In  that  study,  Lg  coda  Q  is 
low  throughout  the  plateaus,  being  in  the  range  150-300.  It  is  higher  for  much  of 
the  peninsula,  ranging  between  350  and  450.  Results  of  this  study  show  that 
is  higher  in  the  peninsTila  than  in  the  plateaus.  The  average  values  for  the 


47 


uppermost  40  km  are  55  and  150  for  Regions  1  and  3,  respectively.  In  the  Lg  coda 
Q  study,  Q  at  1  Hz  varies  laterally  in  the  Arabian  Peninsula  with  higher  values  ( 
about  450)  in  the  shield  area,  and  lower  values  (about  350)  in  the  platform  and 
basins.  In  this  study,  we  obtained  a  model  of  180  for  the  uppermost  20  km  and 
140  for  the  layer  20-80  km  for  the  northern,  western,  and  southern  portions  of 
the  Arabian  Peninsula.  If  the  paths  cross  the  shield,  the  shelf,  and  basins,  we 
have  to  include  a  low-Q  layer  in  the  model  to  explain  the  data.  This  is  also  in 
agreement  with  the  results  from  the  Lg  coda  Q  study. 

9.  Discussion 

Previous  studies  indicate  that  varies  regionally,  being  t3rpically  high  in 
stable  regions  and  low  in  tectonically  active  regions  (Mitchell  1995).  The  northern 
portion  of  the  Middle  East  is  characterized  by  severe  deformation  and  yocmg  tec¬ 
tonic  features.  Continental  collision  is  occurring  along  the  Zagros  suture  zone, 
resulting  in  the  Zagros  folded  belt.  The  left-lateral  East  Anatolian  fault  and 
right-lateral  North  Anatolian  fault  became  active  5  Ma.  As  proposed  in  recent 
research  (Mitchell,  1995;  Mitchell  et  al.,  1997),  and  Q  of  any  wave  propagat¬ 
ing  through  the  crust  in  any  region  are  ultimately  governed  by  the  type,  severity, 
and  temporal  occmrence  of  tectonic  activity  there.  According  to  that  proposal, 
enhanced  temperature  associated  with  tectonic  activity  has  released  water  of 
hydration  that  resides  in  cracks  and  permeable  rock  of  the  crust,  causing  low  Q 
there. 

Tectonically  active  regions  are  usually  characterized  by  high  heat  flow.  Ilk- 
isik  (1995)  recently  found  that  the  average  heat  flow  value  in  western  Anatoha  is 
about  107  mW  m  which  is  about  60%  above  the  world  average  value.  A  value 
as  high  as  247  mW  m~^  was  even  measured  at  a  location  of  38°57'N  and 
29°13.2E.  Wyllie  (1988)  described  ways  in  which  hydrothermal  fluids  could  be 
generated  by  elevated  temperatures  through  broad  regions  of  the  mantle.  In  a 
region  of  plate  convergence,  such  as  Turkey  and  Iran,  those  fluids  may  originate 
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in  a  portion  of  the  Arabian  plate  subducting  beneath  Eurasian  plate.  Upper  man¬ 
tle  fluids,  once  released,  propagate  slowly  upward  and  eventually  reside,  predom¬ 
inantly,  in  the  crust.  The  energy  used  to  move  fluids  through  that  permeable 
crust  is  lost  from  wave  propagating  through  it,  resulting  in  higher  wave  attenua¬ 
tion. 

Although  the  values  in  the  Arabian  Shield  are  higher  than  those  in  the 
northern  portion  of  the  Middle  East,  they  are  lower  than  those  found  for  other 
stable  regions  of  the  world.  Throughout  the  entire  Middle  East  lower  crust 
exhibits  values  no  greater  than  about  150.  This  contrasts  sharply  with  results 
from  other  continental  regions  where  increases  rapidly  with  depth  at  mid¬ 
crust  depths,  usually  having  values  as  high  as  several  hundred  to  a  thousand. 
The  Arabian  Peninsula  has  been  the  site  of  large-scale  uplift  and  extension  with 
associated  volcanic  activity.  The  Afro-Arabian  dome  began  forming  in  the  early 
Cretaceous  and  is  the  result  of  two  or  more  generating  systems  that  are  still 
active  (Almond,  1986).  Temperatures  beneath  the  dome  are  as  high  as  900"C  at  a 
depth  of  40  km  and  these  abnormally  high  temperatures  extend  as  far  as  200  km 
inland  from  the  Red  Sea  (McGuire  and  Bohannon,  1989).  The  low  crustal  Q  in  the 
Arabian  Peninsula  may  therefore  be  due  to  water  released  from  hydrothermal 
reactions  at  those  elevated  temperatures  produced  by  an  upper  mantle  heat 
source. 

10.  Conclusions 

Shear-wave  velocity  models  were  obtained  using  a  two-station  method  for 
three  regions  of  the  Middle  East:  the  Turkish  and  Iranian  Plateaus  (Region  1),  a 
region  surrounding  and  including  the  Black  and  Caspian  Seas  (Region  2),  and  the 
Arabian  Peninsula  (Region  3).  The  shear-velocity  models  of  the  three  regions  dif¬ 
fer  mostly  in  the  uppermost  20  km,  velocity  being  fastest  in  Region  3,  and  slowest 
in  Region  2.  The  velocity  models  for  Regions  1  and  3  agree  well  with  those  found 
in  other  studies  in  the  Middle  East.  Analysis  of  the  resolving  kernels  and  the 
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effect  of  initial  values  on  inversion  indicate  that  our  velocity  models  are  valid  for 
at  least  the  upper  120  km  of  the  crust  and  upper  mantle. 

models  obtained  using  the  two-station  method  show  that  varies  both 
with  depth  and  regionally.  is  higher  in  the  relatively  stable  Arabian  Peninsula 
than  in  the  tectonically  active  plateau  regions  in  the  north  by  a  factor  of  2-3. 
in  the  Turkish  and  Iranian  Plateaus  is  57  in  the  upper  crust,  54  in  the  lower 
crust,  and  89  in  the  upper  mantle,  whereas  values  in  the  Arabian  Peninsula 
at  the  same  depths  are,  respectively,  178,  135,  and  148.  Analysis  of  the  resolving 
kernels  and  the  effect  of  initial  values  on  inversion  indicate  that  our  Qf^  models 
are  only  correct  for  the  uppermost  80  km  of  the  crust  and  upper  mantle. 

Single-station  studies  have  allowed  us  to  study  regional  variations  of 
anelastic  properties  in  greater  detail  than  that  obtained  by  two-station  studies. 
For  Regions  1  and  3,  they  support  our  models  obtained  using  the  two-station 
method.  They  also  show  that  regional  variations  in  occur  within  both  Region  1 
and  Region  3.  The  differences  occur  in  the  uppermost  crust  where  they  vary 
between  50  and  70  in  the  upper  10  km,  between  50  and  60  in  the  depth  range 
10-40  km  for  the  Turkish  and  Iranian  Plateaus,  and  between  140  and  180  in  the 
upper  20  km  of  the  Arabian  Peninsula.  These  results  support  the  idea  (Mitchell, 
1995)  that  in  the  crust  of  any  region  is  directly  proportional  to  the  length  of 
time  that  has  elapsed  since  the  most  recent  major  tectonic  activity  there.  Our 
results  indicate  that  zones  of  plate  convergence  and  regional  upper  mantle  heat¬ 
ing  are  characterized  by  low  extending  from  the  surface  to  upper  mantle 
depths.  This  contrasts  with  results  in  stable  regions  and  regions  of  extension 
where  increases  rapidly  at  mid-crust  depths. 

Average  values  for  both  the  upper  and  lower  crust  in  the  three  regions 
are  lower  than  those  obtained  in  other  continental  regions  of  the  world.  Since  the 
Middle  East  is  a  region  of  plate  convergence  in  the  north  and  high  upper  mantle 
temperature  and  uplift  in  Saudi  Arabia,  fluid-filled  cracks  may  occur  at  lower 
crustal  depths  because  metamorphic  fluids  were  recently  generated  there.  With 
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time,  we  expect  that  those  fluids  will  be  forced  to  shallower  depths  and  at 
lower  crustal  depths  will  become  larger. 
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